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A PROCEDURE FOR COMPUTING 


THE MOTION OF A LUNAR- LANDING VEHICLE 
DURING THE LANDING IMPACT 

By William C. Walton, Jr., and Barbara J. Durling 
Langley Research Center 

SUMMARY 

This paper is a description of a procedure for computing the general motions 
during impact of a vehicle performing a controlled soft landing on the moon. The 
assumptions involved in idealization of the vehicle and landing surface are stated. The 
equations of motion of the idealized system are derived. A procedure is set forth for 
numerical integration of the equations. Detailed and explicit instructions are given for 
coding a digital computer to carry out the computation of an impact history. 

The procedure takes into account experience obtained from extensive correlations 
of computed landing motions with results of tests of a model incorporating many of the 
mechanical actions of a realistic vehicle. These correlations are reported in refer- 
ence 1 (NASA TN D-4215). 

The vehicle is idealized as an arbitrary rigid body with up to four legs arbitrarily 
attached. A leg is composed of three struts in an inverted tripod arrangement. The 
struts telescope during impact. Telescoping of a strut is resisted by forces which 
represent the effects of an aluminum honeycomb shock absorber mounted in the strut 
and the effects of overall system elasticity and damping. Representation of a very gen- 
eral landing surface is allowed. 


INTRODUCTION 

American designers' conceptions of spacecraft capable of controlled, soft landings 
on the moon have been characterized by a preference for vehicles with legs. Prominent 
examples are the Surveyor vehicle (with three legs) for unmanned landing and the Apollo 
lunar module (with four legs) for manned landing. The legs must be designed to prevent 
damage to the craft during the landing impact and to bring it to rest in an upright atti- 
tude so that no part of its mission such as deployment of instruments or relaunch will be 
inhibited. Leg designs have very generally been trusses terminating on the lower end 
in foot pad structures. Load- limiting shock absorbers are built into the trusses. 



Evaluation of the design of such a leg system is a considerable task. The main 
reasons are as follows: 

(1) Limitations on present methods of maneuvering vehicles during the approach 
to landing make it necessary to consider a wide range of attitudes and velocities of the 
vehicle at the moment of impact 

(2) Because of uncertainty regarding the nature of landing sites on the moon, it is 
necessary to consider very general conditions of topography and composition of the 
landing surface 

(3) The relation between the initial conditions and the outcome of an impact gen- 
erally is extremely complicated 

(4) Impact testing of structurally realistic full-scale vehicles in the earth's gravi- 
tational field is seriously hampered by problems which arise in simulating lunar gravity 

(5) Because of fabrication difficulties, it does not appear to be feasible to construct 
models structurally scaled so that the elastic behavior of the model when impacting 
under earth's gravity simulates that of a realistic full-scale vehicle impacting under 
lunar gravity. 

Faced with these difficulties and yet required to make judgements regarding the 
merit of different leg designs, a number of organizations (refs. 2 to 7) have settled on 
the following procedure: 

Equations of motion of an impacting vehicle based on simplifying idealizations of 
the structure and the landing surface are derived. A digital computer program is 
devised for generating numerical solutions of the equations. A dynamic model suitable 
for landing tests on earth is constructed, as many features of a realistic vehicle as is 
feasible being scaled. Landing tests are conducted to determine the behavior of the 
model during impacts under earth's gravity. Impact motion histories are calculated for 
the model and are compared with the results of the tests. On the basis of the correla- 
tion, the analysis is refined; and when consistent success is achieved in predicting the 
behavior of the model, the analysis is used to assess the performance of realistic full- 
scale vehicles impacting under lunar conditions. 

Heavy reliance is currently placed upon such procedures in making design deci- 
sions affecting manned craft. 

The purpose of this paper is to describe in some detail an analytical procedure 
developed at the Langley Research Center along the lines which have been described. 

The model which was used, the testing method, and correlations between the analysis 
and experiment are discussed in reference 1. It is felt that the present paper and ref- 
erence 1 together will be useful as an aid to ascertaining the current state of the art in 
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analytical prediction of the performance of leg systems, in attaining computing capability 
on a level with the state of the art, and in planning research to improve the state of the 
art. 
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SYMBOLS 

time 

space-fixed Cartesian axes 
body-fixed Cartesian axes 
gravitational constant 
Eulerian angles 

space system components of general vector 

body system components of general vector Qx>Qy^Z 

element of matrix transforming space coordinates to body coordinates 
(see eqs. (2) to (3i)) 

body system components of a general vector fixed in body 
space system components of vector C^,C^,C^ 
body system components of angular velocity vector of body 
total mass of vehicle 

principal moments of inertia about body axes for vehicle in its initial 
undeformed configuration 

space system components of center of gravity of body 
body system coordinates of an arbitrary point fixed in body 
space system coordinates of point ijp,»7p,Cp 
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F x ,F y ,F z space system components of total force exerted on body through leg struts 

body system components of total torque about center of gravity of body 
produced by forces through leg struts 

v OX’^OY’ V OZ space system components of velocity of center of gravity of body 

e k kinetic energy of body 

Xp j>Y p j,Zp j space coordinates of foot 

j k’^H j k’^H j k body coordinates of hard point 
•^HJjk’^H j k’^H j k s P ace coordinates of hard point 
instantaneous length of a strut 

^X,j,k’^Y j k’^Z j k space components of unit vector directed along a strut with sense 

from foot to hard point 

F HX,j,k’ F HY j k’^HZ j k space components of force exerted on body through a strut 

F H,j,k magnitude of force exerted on body through a strut, positive when sense is 
from foot to hard point 

Sfi j k amount of contraction of a strut beyond which permanent shortening occurs 

£>R,j,k length of strut when shock absorber is neither contracted nor extended, 

permanent shortening due to crushing of shock-absorbing material being 
taken into account 

S G j ^ initial length of a strut 

Ss j k stroke (see eq. ( 18 )) 

F r j ^ part of magnitude of shock-absorber force termed rate-dependent force 

F<5 • w part of magnitude of shock-absorber force termed quasi-static force 

J? 1 ' 


4 



constants determining rate-dependent force 
(see fig. 2) 

P ; • k coefficient in equation for quasi-static force (see eq. (22)) 

Kg coefficient in equation for stop force (see eq. (23)) 


S ECJ,k> s EE,j,k’ S CCJ,kScE,j,k 


F EC,j,k’ F EE,j,k’ F CC,j,k’ F CE,j,k 


F H£,j,k’ F HT?,j,k’ F HC,j,k body com P° nents of force F HX ^ k ,F H Yj >k ,F HZ j jk 


L X,j’ A Y,j’ A Z,j’*j 


coefficients in equation describing landing- surface boundary 
plane associated with jth foot (see eq. (27)) 


W v -,W 7 normalized coefficients in equation describing landing- 

surface boundary plane (see eqs. (28) to (29c)) 


Hj signed distance of a point from jth landing surface plane (see eq. (30)) 

Np qj,Tp q j elements of matrices for computing components of a vector normal 

and tangential to a boundary plane, respectively (see eqs. (31a) 
to (32b)) 


F FNSX j ,F FNSY,j ,F FNSZ,j P ar t °f force upon a foot acting normal to landing surface 

termed quasi- static normal force 

FpNS j magnitude of quasi- static normal force 

Kj j coefficient in equation for Fpjqg j (see eq. (38)) 

Dj absolute distance from foot to landing surface plane (see eq. (39)) 

F FNDX j ,F FNDY j’ F FNDZ j P ard °* f° rce u P on a foot acting normal to landing surface 

termed dynamic normal force 

Rj^ j coefficient in equation for dynamic normal force (see eq. (40)) 

XpN i’^FN i’^FN i space components of velocity of foot normal to landing 

,J ,J surface plane 
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F FTDX,j F FTDY,j’ F FTDZ,j force on a foot acting tangential to landing surface plane 

termed dynamic tangential force 

R T j coefficient in equation for dynamic tangential force (see eq. (42)) 

^FT, j’^FT, j’^FT, j space components of velocity of foot tangential to landing 

surface plane 

F FLX,j ,F FLY j’ F FLZ j space components of total force on a foot acting through 

three struts bearing on foot 

F FLNX j ,F FLNY j ,F FLNZ j s P ace components of component of force 

F FLX,j’ F FLY,j F FLZ,j normal to landing surface 
plane 


F FLTX,j> F FLTY,j’ F FLTZ,j space components of component of force 

F FLX,j’ F FLY,j’ F FLZ,j tangential to landing surface 
plane 


t (f) 


some time at which a foot has just become free of landing surface 


-yV-p ; body coordinates of foot 

” »J * .1 ” , 1 


t(f) n (D r (f) 
W,J F,j 

X^ ,Y^ ,Z® 
F >1 F,j F,j 


body coordinates of foot becoming free at time t^ 
space coordinates of foot becoming free at time t^ 


L some characteristic length of the vehicle 


The subscript B on any symbol denotes a nondimensional value. A dot over a 
symbol denotes differentiation with respect to time t. An asterisk over a symbol 
denotes differentiation with respect to tg. The superscript (n) denotes nth iteration. 
| | denotes absolute values. 
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GENERAL ASPECTS OF IDEALIZATION 


Overall Vehicle 

The vehicle is treated as an arbitrary rigid body to which there are attached up to 
four legs, each leg consisting of three struts in an inverted tripod arrangement. (See 
fig. 1.) The struts are connected to the body by universal joints, and the junction point 
of the three struts in a leg is also a universal joint. This junction point is called the 
foot of the leg, and the points where struts attach to the body are called hard points. 
There is a shock absorber in each strut. The individual struts may shorten or lengthen 
because of stroking of the shock absorbers, but otherwise they do not deform. Locations 
of the hard points on the body and the initial positions of the feet relative to the body may 
be arbitrarily chosen. 



Figure L- Idealization of vehicle. 

Inertias 

In the idealized system the legs are considered to have no mass. The inertial 
properties of the body are characterized by specification of the total mass, the center of 
gravity, a set of principal axes, and the moments of inertia taken about the principal 
axes. In representing an actual vehicle or model, the inertial properties are computed 
for the system as a whole, including the legs. These inertial properties are then 
assigned to the body alone in the idealized system. This approach has been adopted 
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because it permits considerable simplification of the equations of motion and reduces 
the number of equations. The approach is reasonable if the legs are of light construc- 
tion and/or deformations of the legs due to stroking of the shock absorbers are small; 
that is, the moments of inertia do not change appreciably because of stroking of the 
struts. 


Shock Absorbers 

The shock absorber in a strut is considered to produce simultaneously a force at 
the hard point and at the foot to which the strut is connected. The forces are considered 
to be equal in magnitude but opposite in direction and to be directed along the axis of the 
strut. The magnitude of the shock absorber force is assumed to depend on the change 
of length of the strut from the initial length and on the rate of change of the length. The 
specific relation between the force and these variables is described in the section 
"Forces and Torques on the Body.” 


Landing Surface 

The boundary of the landing surface is represented by a set of arbitrarily oriented 
planes, one plane associated with each foot. Use of a different plane for each foot allows 
for the representation of many irregular surfaces. When the foot of a leg is not inter- 
acting with the landing surface material, the leg is assumed to move as a rigid extension 
of the body. When the foot is interacting with the landing surface material, a force acts 
on the foot, and the shock absorbers may stroke. The force exerted on the foot by the 
landing surface is assumed to be a function of the position and velocity of the foot. A 
basic assumption of the analysis, leading to equations of motion of the feet, is that a foot 
always moves in such a way that this force is kept in balance with the shock-absorber 
forces bearing on the foot. 

In representing the forces generated by the interaction of the foot pad with the 
landing surface material, the analyst must face two facts. First, knowledge of the 
properties of the lunar surface which would affect landing performance is as yet limited. 
Second, soil mechanics has not progressed to the point where one can predict with any 
confidence the history of forces on an arbitrary body impinging upon or passing through 
soil even under laboratory conditions. The reaction to these difficulties has been a 
rather general concentration of initial analytical effort on the case where a foot upon 
contacting the surface is stopped abruptly, is effectively pinned, and remains in place 
until there is a tendency for it to lift off. Cases where the feet move substantially 
through or along the surface are studied by assuming laws for the force on a foot. The 
procedure given by this paper is organized so that one may construct a variety of laws 
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for the force on a foot by specialization of constants. The details of the computation 
are given in the section "Equations of Motion of the Feet." 

EQUATIONS OF MOTION OF THE BODY 

The object in this part of the paper is to set down the differential equations which 
govern the motion of the central rigid body. Equations of motion of the feet are devel- 
oped in a subsequent section. 

Coordinate Systems 

Space and body coordinate systems .- Reference is made to space and body coor- 
dinate systems (fig. 2). These systems are both right-handed Cartesian systems. The 
space system is assumed to be an inertial 
system, and its axes are denoted by X, Y, Z 

and Z. It is oriented so that gravitational 
forces point along the negative Z-axis. The 
body system is assumed to be fixed in the 
body, and its axes are denoted by £, r 
and It is placed within the body so that 
the origin coincides with the center of grav- 
ity and the three coordinate axes coincide 
with principal axes of inertia. Further con- 
siderations explained under the heading 
"Eulerian Angles" also enter into the estab- 
lishment of the body axes. 

Eulerian angles.- Rotations of the body X 
are specified by use of Eulerian angles 0, Flgure 2 _ Space and ^ coordinate system 

6, and 0/ as shown in figure 3. (The hori- 
zontal orientation of the Z-axis in the sketch is merely for convenience in drawing.) 

Note that the Eulerian angles are not defined in the usual way; that is, setting the angles 
equal to zero does not bring the body axes into coincidence with the space axes. Instead, 
the positive £-axis merges with the positive X-axis, the positive ? 7 -axis with the negative 
Z-axis, and the positive £-axis with the positive Y-axis. 

The following limits are set on the ranges of the Eulerian angles: 


0 < 6 < 7T 

(la) 

VII 

•e- 

VII 

(lb) 

VII 

VII 

(lc) 
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Figure 3.- Eulerian angles. 


The limits on <p and xp do not restrict 
the generality of orientations. The limits 
on 0 are restrictive, however, and are 
imposed to avoid encountering a mathemat- 
ical singularity for 6 equal to 0 or u as 
explained in a subsequent section. As 6 
approaches either of the end values 0 or tt, 
the t^-axis approaches a horizontal orienta- 
tion. Therefore, in positioning the body 
axes within the body, the £-axis should be 
associated with the principal axis of the 
vehicle which is most nearly longitudinal. 

If this is done, the restriction on 6 causes 
no practical limitation on motions, because 


calculation normally stops before the vehicle has tipped so far that the longitudinal axis 
becomes horizontal. 

Transforming vector components between space and body systems .- Let Q-^, Qy, 
and be the space system components of a vector, and let Q^, Q r ^ ? and be 

the body system components of the same vector. As described in reference 8 and other 
texts on classical mechanics, the components are related by the equation 


N 






s, - 

i 


6 p,q 

< 


L 

-j 

i 

i \ 


(p = 1, 2, 3; q = 1, 2, 3) 


( 2 ) 


Qz 




where 


H - cos xp cos (p - cos 6 sin (b sin xp 

(3a) 

= sin xp sin 6 

(3b) 

13 = “ cos ^ s ^ n “ cos 6 cos ( P sin ^ 

(3c) 

2 i = -sin xp cos < p - cos 0 sin <p cos xp 

(3d) 

22 = cos xp sin 6 

(3e) 

'23 = sin xp sin (p - cos 6 cos 6 cos xp 

(3f) 

>31 = sin 6 sin (p 

(3g) 

'32 ” cos 0 

(3h) 

'33 = sin 6 cos (p 

(3i) 
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The matrix jS p ^ appearing in equation (2) is orthogonal which means that its inverse 
and its transpose are identical. Therefore, the inverse transformation giving space 
components in terms of body components is 


r 

Q x 



T 

r 

<Qy 

s .... 

/ 

6 p,q 


Q?7 ) 

Qzj 




Qr 

L V 


the T denoting the transpose. It follows from equation (4) that if C^, C^, and 
are the constant body system components of a vector fixed in the moving body, and 
Cy, an< 3 Cg are the space components of the vector then 


(4) 


C X> 


r. ~s 

C X 



T 

r 

< 

) = 

^p,q 

< 

c i i ) 


where a dot denotes differentiation with respect to time. 

Relation between Eulerian angles and direction cosines .- The elements of the 
matrix |6 p ^qj in equation (2) may be interpreted as the direction cosines of the body 
axes measured with respect to the space axes. That is, if one imagines a translation of 
the space axes without rotation such that the origins of the body system and the translated 
space system are brought into coincidence, 6^ is the cosine of the angle between the 
£-axis and the translated X-axis, 6^ is the angle between the £-axis and the trans- 
lated Y-axis, and so forth. 

It is generally considerably easier to specify the direction cosines associated with 
an initial orientation of the body than it is to specify the Eulerian angles. Therefore, 
it is useful to have a procedure for computing Eulerian angles once direction cosines 
are given. The following scheme which is easily derived from equations (3a) to (3i) has 
proved to be satisfactory. 


°33 

|$| = arc cosine -■■■ 

+ P33) 2 

VII 

]e[ 

VII 

o. 

(6a) 

II 

(631 g 0) 

(6b) 

1 

II 

-e- 

(631 < 0) 

(6c) 

6 = arc cosine 632 

(0 < 6 < tt) 

(6d) 
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i i 60? 

|0| = arc cosine — 

Vil 

3 E 

VII 

3 

(6e) 

J/c \2 /- \2 

31 ) + ( 6 33 ) 



II 

IT 

(«12-°) 

(6f) 

ii 

( 5 12 < °) 

(6g) 


Relation between Eulerian angles and components o f angular velocity.- Let 
oo and u> ^ denote the body components of the angular velocity vector of the body. 
They are related to the rates of change of the Eulerian angles as follows: (See p. 134 

of ref. 8.) 


r “a 


- *1 

r. 'n 



sin 0 sin xp cos xp 0 1 

i 0 


- 

sin 0 cos xp -sin xp- 0 

In 

00^ 


cos 0 0 1 

0 


The inverse relation is 


f .'s 


sin xp 

COS 0 

0 



1 


sin 0 

sin 0 




> = 

cos xp 

-sin 0 

0 

i 


ip 

1 

j 

cos 6 sin xp 

COS 9 COS 0 

1 


c Or 

V ) 


_ sin 0 

sin 0 

• 


s 

v. J 


The singularity referred to previously in imposing a restricted range on 0 is 
involved in the factor 1/sin 0 appearing in equation (8). With the range as established, 
the denominator sin 0 never vanishes. 


Equations of Motion 

The following definitions are necessary for expressing the equations of motion of 
the body: 

M total mass of vehicle 


x o> Y o> z o 

£p ,7 7p f ^P ''l 
x p, Y p,Zpj 


principal moments of inertia about body axes £, ? 7 , and £ 

space system components of center of gravity of body 

body and space system coordinates, respectively, of an arbitrary point fixed 
in body 
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F x ,F y ,F z space system components of total force exerted on body through leg struts 


body system components of total torque about center of gravity of body 
produced by forces through leg struts 


^OX» v OY>^OZ space system components of velocity of center of gravity of body 

g gravitational constant 

By using these definitions, the elementary equations of translation are 


x o = Vox 


(9a) 

Vo = V OY 


(9b) 

= V oz 


(9c) 

M V 0 X - F x 


(9d) 

MV o Y ~ 


(9e) 

MV oz = F z - Mg 


(9f) 

Euler’s equations for rotation of the body are 



= ( w t?)( w ?)( l7 » ' 

) +N{ 

(10a) 

= ( a, e)( a, s)( I c _ h) 

| + 

(10b) 

= ("ONX 1 * ■ H 

+ 

(10c) 


For completeness, equation (8) is repeated 

r ~\ r 




<p 

*> = 


V J 


sin \f/ 
sin 0 

cos i p 


cos j// 
sin 6 

-sin \fs 


cos 8 sin j// cos 8 cos \p ^ 
sin 8 sin 8 


r ^ 
c o 


CjO 


ri 


CO ^ 

v y 


(ii) 
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From equation (4), it follows that 


X P - Xq 



H 

S' 

V 

i 

>< 

o 


6 p,q 


Vp 

Zp - Zq 

V V 

! 



Cp 

J 


( 12 ) 


The preceding equations constitute the equations of motion of the body. A dimen- 
sionless form of the equations more suitable for use in computation is given in a sub- 
sequent section. 


Kinetic Energy 

It is sometimes useful to compute a history of the kinetic energy of the system 
throughout an impact. The kinetic energy is given by the equation: 

% = f( v ox + v oy + v o^ + + 1 n U) v + I ? w ?) 

FORCES AND TORQUES ON THE BODY 


The forces which act upon the central rigid body are the gravitational force and 
the forces through the leg struts produced by the action of the shock absorbers. The 
object here is to write down rules for computing these forces based on practical con- 
siderations which arose in analyzing the behavior of the model of reference 1. 


Definitions 

It is assumed that the legs have been numbered 1 to 4 in any order and the struts 
in each leg have been numbered 1 to 3 in any order. Cases where there are less than 
four legs are handled by setting the shock absorber forces equal to zero in some legs. 
The following definitions are necessary to subsequent developments: 


X F,j’ Y F,j’ Z F,j 


instantaneous space coordinates of foot of jth leg 


X H,j,k’ Y H,j,k’ Z H,j,kJ 


body and space coordinates, respectively, of hard point at which 
the kth strut of jth leg is attached to body 


S. , instantaneous length of kth strut of jth leg 

J,K 


U Y i i,,U v . . , U rr • i space components of a unit vector directed along kth strut of 

jth leg, the sense being from foot to hard point 
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F HX,j,k’ F HY,j,k’ F HZ,j,k 


space system components of force exerted on body 
through kth strut of jth leg 


F H,j,k magnitude of force F HX,j,k’ F HY,j,k’ F HZ j k taken positive if the force is 

directed from foot to hard point and negative if the force is directed from 
hard point to foot 

The following equations may be written on the basis of these definitions: 


S j,k ~ _( X H,j,k - X F,jf + ( Y H,j,k “ Y F,j) 2 + ( Z H,j,k ‘ Z F,j)' 


1/2 


( 14 ) 


3j ’ k \k ( Xli ’j’ k ’ X F,j)( X H,j,k ' X F,j) + ( Y H,j,k “ Y F,j)( Y H,j,k ' Y F,; 
+ ( Z H,j,k- Z F,j)( Z H,j,k- Z F,j) 


( 15 ) 


U XJ,k’ U Y,j,k’ U Z,j,k-s^ 


X H,j,k “ X F ,i)> ( Y H,j,k “ Y F,j)» ( Z H,j,k “ Z F,j 


( 16 ) 


F HXJ,k’ F HY,j,k’ F HZ,j,k 


= F H,j,k( U X,j,k ,U Y,j,k ,U Z,j,k) 


( 17 ) 


Reference Length 

The shock absorbers for the model of reference 1 were cylinders of aluminum 
honeycomb mounted in the struts. For compressive loading this type of shock absorber 
has a load-deflection relationship as shown in figure 4. 
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The shock absorber compresses because of elastic deformation until the load 
reaches the level at which crushing of the honeycomb begins. Compression continues at 
constant load as the cylinder crushes. The elastic part of the load-deflection diagram 
may be nonlinear. If crushing occurs, the cylinder becomes permanently shortened, 
and leaves a gap as shown by figure 5. Actually the shape of the elastic part of the load- 
deflection diagram may be altered somewhat as a result of crushing; but in this analysis, 
it is assumed that this shape remains unchanged. 



In order to account for permanent shortening, it is convenient to define a variable 
called the reference length of a strut denoted by the symbol Sp • Initially, the ref- 
erence length is equal to Sq j ^ the initial length of the strut. The stroke Sg j ^ at 
any time t is defined to be the reference length minus the length of the strut at time t. 

S S,j,k = S R,j,k " S j,k < 18) 

The length of the portion of the deflection denoted ’’elastic'' in sketches 3 and 4 is repre- 
sented by the symbol S E - At the completion of any stroke exceeding S E . the 
reference length is diminished by the amount by which the stroke exceeded S E j 


Magnitude of Shock Absorber Force 

It is convenient to express Ftt i the magnitude of the shock absorber force, as 
the sum of two terms 


F H,j,k F R,j,k + F S,j,k 


(19) 


The term F D . , is spoken of as the rate-dependent force and F c . , as the quasi- 

b )J> K 

static force. 


Rate-dependent force .- Figure 6 shows the relation between the rate -dependent 
force and the dynamic state of the strut. With respect to the current reference length, 
a strut is either extended j k - Sj ^ = Cfj or contracted j k " 8 j k > 0). 1° either 
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(a) Extended strut, Srj^ - Sj^ = 0. (b> Contracted strut, Srj |< - Sj ^ = 0. 

Figure 6.- Relationship between rate- dependent force and the dynamic state of the strut. 

case the strut may be extending (Sj >k > o) or contracting (Sj^ ^ oj. For each of the four 
possible combinations of these conditions, the rate-dependent force F R . k is described 
by a separate ramp function of Sj jk , the rate at which the strut length is changing. Eight 
arbitrary constants associated with each strut establish the shapes of the ramps. These 
constants are S E q j k > ^EE j k J ^CC j k* a ^ d ^CE j k which bear the units of velocity 

and F ec j >k , F EE,j,k’ F CC,j,k> and F CE,j,k which bear the units of force. As fig- 
ure 6 indicates, the velocity constants establish ranges of Sj jk within which F R j k 
varies linearly with Sj >k and beyond which F R j k is a constant. The force constants 
establish the heights of the ramps. The subscripts EC, EE, CC, and CE, respec- 
tively, denote the conditions extended and contracting, extended and extending, con- 
tracted and contracting, and contracted and extending. 

Quasi-static force. - The quasi-static force is considered to be zero unless one of 
the following conditions holds: (1) The strut is contracted with respect to the current 

reference length and is extending 


! R,j,k " S j,k > 0 

(20a) 

! j,k > 0 

(20b) 
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(2) The strut is extended beyond its initial length S 0,j,k and is extending 


S j,k - s 0,j,k = 0 


S j,k > o 


(21a) 

(21b) 


When the first condition holds, the force is assumed to be given by the following 
third-order polynomial function of the stroke Sp^ j k " k : 

F S,j,k = P l,j,k + P 2,j,k( S R,j,k - S j,k) + P 3,j,k( S R,j,k - S j,k) 2 + P 4,j,k ( S R, j,k " S j,k) 

( 22 ) 

wherein the constants P 7 . , (l = 1, 2, 3, 4) may be arbitrarily selected. For this con- 
dition the stroke is always less than or equal to Sg j k because of the way S R ^ k is 
reset at the completion of a period of contracting which necessarily precedes a period 
of extending. Figure 7 shows an example of how the quasi-static force might vary with 
the stroke S R>j k - S j k . 



S o , - S. , 
R , j , k j , k 


Figure /.- Quasi-static force diagram for strut contracted but extending. 


When the second condition holds, the equation for the quasi-static force is: 

F S,j,k = " K s( S j,k - S OJ,k) (23) 

where Kg is an arbitrary positive constant. 

Representing Forces for the Model of Reference 1 

The preceding relations between the magnitude of the shock absorber force and the 
dynamic state of a strut were established with the following three objectives in mind: 
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(1) to represent constant force crushing and free extension characteristic of the alumi- 
num honeycomb shock absorbers used on the model of reference 1, (2) to represent 
measured elastic characteristics of the model, and (3) to represent, in a gross way, 
energy dissipation due to causes other than crushing of the honeycomb to prevent unreal- 
istic bouncing of the idealized system. 

Constant force crushing .- To accomplish the first objective, the constant F cc j k 
of the rate-dependent force is set equal to the force at which the honeycomb crushes, 
and the constants ^EE,j,k’ an< ^ ^CEAk are ass ig ne d relatively small values 

for a nominal representation of frictional resistance to stroking produced by the bearings 
in the strut. The constants Sgcj )k , Sgg j k , Sqqj k , and are assigned 

values greater than zero but small compared with the average stroking rate of the shock 
absorbers. With these settings, the shock absorbers will operate on the constant-force 
parts of the rate-dependent force diagrams except when stroking rates are very low. 

When stroking rates are very low, it is desirable to operate on the sloping part of the 
rate-dependent curves in order to avoid an annoying problem in numerical integration of 
the equations of motion. The problem is an oscillatory instability of the calculated 
stroking rates Sj k which is encountered when very large forces occur simultaneously 
with very small stroking rates so that the sign of a stroking rate may change on one time 
step and change back on the next time step. 

Elasticity .- For the second objective, the constants Sg j k and j k are 
chosen for all struts to represent insofar as possible the elastic behavior of the model. 
Since there can be significant elastic deformations of the body as well as of the struts 
and shock absorbers, the validity of such a representation is open to question. This is 
an important limitation arising out of the simplifications made in the analysis. Refer- 
ence 1 describes in detail a procedure used to assign these constants for the test model. 
The procedure involved a series of static load tests of the model with various leg struts 
removed. In spite of the questionable aspects of the representation, the procedure 
resulted in good predictions of the outcome of test impacts in which elastic rebound was 
a significant influence. 

Damping .- It is noted that for the idealized system, there may be an axial elastic 
force in a strut which is contracted with respect to the reference length only when the 
strut is extending. If a strut is contracted with respect to the reference length and is 
contracting the full crush force, F cc • k resists the stroke unless the stroking rate is 
less than the constant S C( ~. j k . As discussed previously, Sqq j k is set very low com- 
pared with the stroking rates expected throughout most of the impact. The elastic force 
must be lower than the crush force. Therefore, for any appreciable loading rate, it 
takes more energy to compress a strut spring than is given up when the spring unloads. 

As a result, elastic action is damped. This representation is not realistic. A strut 
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containing a crushable aluminum honeycomb shock absorber both loads and unloads 
elastically with a period of crushing between if the loads get high enough. Energy losses 
other than those caused by crushing of the honeycomb come about in diverse ways. How- 
ever, the representation described here has been found to be an easy way to provide 
some damping for the idealized system during that part of an impact history during which 
loads have subsided to the point where the honeycomb in the shock absorbers no longer 
crushes. If such damping is not provided, the idealized system will bounce erratically 
throughout this period. This result is contrary to the behavior of the model of refer- 
ence 1 during test impacts. There the motion was quickly damped to that of a smooth, 
essentially rigid body, pivoting about one or two feet pressed against the landing surface. 
Unrealistic bouncing can result in erroneous predictions as to whether the model will 
overturn as the result of an impact. Representation of the elasticity of the model is not 
marred by removing the elastic forces from the compressive part of the stroke. Elastic 
action of the model appears to be significant only during a lew initial relatively hard 
impingements of the feet on the landing surface. On the compressive part of the stroke 
during such an impingement, the strut compresses through the elastic part of the stroke 
so quickly that it does not matter what representation of the force is used. The primary 
effect of elastic forces comes as they impart a push to the vehicle during extension of a 
strut. 

Quasi-static force for an extended strut .- The force given by equation (23) is not 
considered to be of very general interest. It was included in the programing in order to 
keep struts from extending beyond their original lengths since the struts of the model of 
reference 1 could not do so. A force which was linear in the extension rather than a 
step force was chosen to avoid instability in the numerical integration of the type dis- 
cussed previously. The force acts only if the stnit is extending, and it opposes exten- 
sion. Consequently, it can only absorb energy and can never add energy to the system. 
This aspect of the programing is important because the constant Kg must be set very 
high to prevent substantial extension of the strut. If the force is programed as a spring 
which can return energy to the system, lightly damped high-frequency oscillations may 
be introduced into the calculated motion and cause difficulties in the numerical 
integration. 


Total Forces and Torques 


The space system components of the total force acting on the body are computed 
by the equation 


j=4 k=3 

F X’ F Y’ F Z = 2 Z ( F HX,j,k’ F HY,j,k F HZ,i,k) 
j=l k=l 


(24) 
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The body system components of the force acting on the body at a hard point are 

denoted by F H| ) j ) k’ F H? 7 ,j,k ,F HC,j,k and are related to the space system components 
of the force by the equation 


r z 


F H?7,j,k 

= 

F H?,j,k y 



HX,j,k 
{ F HY,j,k 
F HZ,j,k 




( 25 ) 


The body system components of the total torque acting on the body are given by 
the equations: 



j=4 

k=3 




II 

.0^ 

1 (^H, j,k F H^, j,k ‘ 

" ^H,j,k F H?7,j,k) 

(26 a) 


1=1 

k=l ' 



j=4 

k=3 



N 

!W 

11 

) ( C H,j,k F H£,j,k ‘ 

‘ ^H,j,k F H?,j,k) 

(26b) 


]=1 

k=l 




& 

k=3 



N* 

~ z 

Z (^H,j,k F Hrj,j,k " 

^H,j,k F H^,j,k) 

(26c) 



k=l 



EQUATIONS OF MOTION OF THE FEET 


Considerations Related to Landing Surface Planes 

Equations of planes .- The landing surface plane associated with the jth foot is 
described by an equation of the form: 


A X,j X + A Y,j Y + A Z,j Z+ V° ( 27 ) 

where A x ^, A y ., A z and A^ are constants which may be artibrarily selected. 

The equation is rewritten with normalized coefficients as follows in order to avoid 
ambiguity in specifying the position of a point in space relative to the plane: 


W v .X + W v .Y + W„ .Z + W. = 0 
A >J Z,J ] 


(28) 


where 


W X J j’ W Y,j’ W Z,j’ W j = 


A xj 


A X,j’ A Y,j’ A Z,j’T7 


Z, j\| ( A X,j) + ( A Y,j) + ( A Z,j) 


( A z,j * °) < 29a ) 
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Normals.- The vector W v .,W V .,W 7 . is normal to the landing surface plane 



associated with the jth foot. As a result of the normalization procedure just described, 
the normal is a unit vector directed as follows: 

(1) If the jth plane is not vertical, the projection of the normal on the Z-axis points 
in the positive Z -direction. 

(2) If the jth plane is vertical but not parallel to the Y-axis, the projection of the 
normal on the Y-axis points in the negative Y-direction. 

(3) If the jth plane is vertical and parallel to the Y-axis, the normal points in the 
negative X-direction. 

Surface and subsurface sides of a plane .- A landing surface plane divides inertial 
space into two spaces, one into which the normal is directed and one out of which the 
normal is directed. Points lying within the space into which the normal is directed are 
said to be to the surface side of the plane. Points lying within the space out of which the 
normal is directed are said to be to the subsurface side of the plane. 

Distance of a point from a plane .- Let the symbol Hj denote the length of the per- 
pendicular from a general point X,Y,Z to the jth landing surface plane. The distance 
is considered positive if the point is to the surface side of the plane and negative if the 
point is to the subsurface side of the plane. The following formula then gives the 
distance: 

Hj = W X j X + W Y j Y + W Zj] Z ♦ Wj (30) 

Normal and tangen ti al projections of a vector .- If the space system components 
Qx>Qy’ Q Z an ar ^iti'ary vector are given, it will be necessary to resolve the vector 
into the sum of a vector normal to and a vector tangential to the jth landing surface plane. 
The components in the space system of the normal and tangential vectors are, respec- 
tively, N^r j?Ny j,N z j and T x j? T y j> T z j and are readily computed with use of the 
following relations: 
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N X,j 




r "\ 

Q x 

< N Y,j 

> " 

N 

p,q,j 

< 

| QY i 

N Z,j 


- 


Qz 

v J 

r *\ 

T X,j 

i 

| 

r 


c 

Q x 

< T Y,j 

\- 

T 

A p,q,j 

< 

Qy^ 

T z,i 

V J 




Q z 

^ J 


where 



2 

'h 

i 

ii 

( W X,i)( 


( W X,i)( 

] W Y,i) 2 ♦ 

( w z,i ) 2 

-( w x,j)l 

: w y,j) 

-( w x,j)l 

( w z.j) 


Ki)( W Y,i) 
( W Y,i) 2 
Z,j) 

-( W X,i)(WY,i) 

( w x,j) 2 + (w Zji ) 2 

-( W Y,i)(*Z,j) 



( w x,j)( w z,j) 

( W Y,i)( W Z,i) 

( w z,i ) 2 

-(Wx,i)(Wz,i) 

-( W Y,i)( W Z,i) 
|( W X,j ) 2 ♦ (WY,i ) 2 


(31a) 


(31b) 


(32a) 


(32b) 


Free and penetrating foot .- The following definitions are made to facilitate concise 
statements regarding whether a foot is interacting with the landing surface material. 

The jth foot is said to be free if any of the following conditions hold: 

(a) The foot is on the associated landing surface plane 


w X,j x F,j + W Y,j Y F,j + w Z,j z F,j + Wj = 0 

(b) The foot is to the surface side of the plane 

W X,j X F,j + W Y,j Y F,j + W Z,j z F,j + Wj > 0 


(33) 


(34) 


(c) The foot is to the subsurface side of the plane; and at the same time, the veloc- 
ity component of the foot normal to the plane tends to carry the foot toward the plane 
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(35a) 


The jth foot is 
the same time 
tends to carry 


w x,jX F>i ♦ w Y>j y Fjj ♦ w z>i z F>j + Wj < o 


W X,i X F,i +W Y,i Y F,i + W Z,fF,i>» 


(35b) 


said to be penetrating if it is to the subsurface side of the plane, and at 
the velocity component of the foot normal to the plane is either zero or 
the foot away from the plane 


W X,j X F ,i + W Y,i Y F ,j + W Z,j Z F,j + W j<° 


Wx,i X F,i + W Y,j Y F,j * W Z,j Z F.l ? 0 


(36a) 

(36b) 


Forces on a Penetrating Foot 


When a foot is penetrating, it is considered to interact with the landing surface 
material and thus produces a force on the foot. Additional forces act on the foot through 
the three struts bearing on it. The force caused by the interaction with the landing sur- 
face material is considered to be the resultant of three separate forces termed the 
quasi-static normal force, the dynamic normal force, and the dynamic tangential force. 


Quasi-static normal force.- This force acts normal to the landing surface plane 
associated with a foot and may therefore be expressed for the jth foot as 

F FNSX,j’ F FNSY,j’ F FNSZ,j = F FNS, j ( W X,j’ W Y,j’ W Z, j) 


(37) 


Here FjtnsX y ^FNSY j 5 anc * ^FNSZ j are s P ace system components of the force, 

F^xto • is a scalar, and W v W v and . are the components of the unit nor- 

FNS,J ’ x,]’ Y,y Z,j 

mal vector previously defined. The scalar F FN g . is taken to be a cubic function of 

the variable Dj which is defined to be the absolute length of the perpendicular from the 

jth foot to the associated landing surface plane. 


F rH5 . s K. .D. + K 0 .D 2 + Ko -D 2 
FNS,] 1,] ] 2,j j 3,3 ] 


(38) 


where the coefficients K; ; (i = 1, 2, 3) may be arbitrarily assigned. From equa- 

A > J 

tion (30), Dj may be computed by the formula 


D. 


W X,3 X F,3 


+ W Y,3 Y F,j-W Z)j Z F)j + W3 


(39) 


Dynamic nor mal force .- The space system components of this force are denoted 
F FNDX j’ F FNDY j’ anc * f FNDZ j- ^he force is proportional to the component of the 
velocity of the foot normal to the associated landing surface plane and acts in the direc- 
tion opposite to that of the normal velocity; that is, 
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r FNDX,j’ F FNDY,j’ F FNDZ,j = l X FN,j» Y FN,j’ Z FN,j) 


(40) 


where Rjjj is a positive constant to be arbitrarily assigned and X FI ^ j, Yp^ j, and 
ZpN,j are the space system components of the normal velocity of the foot. By 
equation (31a), 




\ Y FN, j 

> = 

, Z FN,jj 



N. 


p,q,i 


X F,i 

*F,J 




(41) 


The dynamic normal force has the character of a frictional resistance to penetra- 
tion into the landing surface by the foot. The constant R N j has the character of the 
reciprocal of a viscosity constant. 

Dynamic tangential force .- This force is the tangential counterpart of the 
dynamic normal force. The space system components of the force, denoted 
F FTDX,j’ F FTDY, j ,F FTDZ,j’ are §^ ven by the equation 


' FTDX,j ,F FTDY,j ,F FTDZ,j “ B^^FT^FT^FT,]) 


(42) 


where R*p j is a positive constant, and Xprp j* *^FT j* aa tl Zpip j are the space sys- 
tem coordinates of the tangential velocity given according to equation (31b) by 


(43) 


X FT,j 

! 



X F,j 

<JF T,j 

II 

T p,q,j 

< 

*F,j> 

^FT,j 




Si 


Force through the struts .- The force acting on a foot through a strut is equal in 
magnitude but opposite in direction to the force F HX j j k ,F HZ j k which is 

exerted through the strut onto the hard point at which the strut is attached to the body. 
Therefore, the space components F FLX j, F FLY j, and F FLZ j of the total force on 
a foot through the three struts bearing on the foot may be computed by the formula 


k=3 


Ftht v i.Frit t • . t ' Tn T n • — / F jxv i i.. Fmr ■ ,_.F T 

J >] Lj 


FLX,j’ FLY,j’ FLZ, 


k=l 


HX,j,k’ HY,j,k’ HZ,j,k 


(44) 
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By using equations (31a) and (31b), this force can be resolved into a force normal to and 
a force tangential to the associated landing surface plane with space components given 
by the respective equations: 


f] p ''l 

FLNX,j 





F FLX,j 

F FLNY,j 

= 


N 

p>q,3 

< 

F FLY,j 

F FLNZ,j 





F FLZ,j 


> 



- 


Iw. 

J 




p 



( U 


f FLTX,j 





FLX,j 

F 

*FLTY,j 



T 

A p,qj 

< 

F FLY,j 

F FLTZ,j 





F FLZ,j 

V J 


(4 5a) 


(45b) 


Equation for a Penetrating Foot 

Summing the forces on a foot in the normal and tangential directions gives the fol- 
lowing equations: 


X FN,j’ Y FN,j’ Z FN,j " F FNS,j R N,i( W X,j ,W Y,j’ W Z,j) + R N,i( F FLNX,j» F FLNY,j»' F FLNZ,j 

(46a) 

(46b) 


X FT,j’ Y FT,j’ Z FT,j - R T,i ( F FLTX,j’ F FLTY,j’ F FLTZ,j j 


Adding equations (46a) and (46b) gives 

X F,j’ Y F,j’ Z F,j = F FNS,j R Nj( W X,j’ W Y,j’ W Z,j) + R N,j ( F FLNX,j’ F FLNY, j’ F FLNZ, 


+ R T,j ( F FLTX,j’ F FLTY,j’ F FLTZ,jj 


(47) 


Equation (47) is the equation of motion for the jth foot when the foot is penetrating. 


Equation for a Free Foot 

For computing the trajectory of a free foot, it is assumed that the lengths of the 
three struts connecting the foot to the body remain fixed so that the foot moves as a 
rigid extension of the body. Let . denote the instantaneous body coordi- 

nates of the jth foot. From equation (2) 
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Suppose that at a specified time D f ' a foot has just become free. This time may be 

the initial time or a later time following a period of penetration. The space coordinates 

of the foot at this time may be presumed to be known and are denoted by ,Y^ 

F,f F,j’ F,j‘ 

The corresponding body coordinates are denoted by 7 /^., and . and may be 

1 F,j F,j 


The corresponding body coordinates are denoted by and . and may be 

F,j F,j 

computed from equation (48). The body coordinates of the foot remain fixed as long as 
the foot remains free. Therefore, during this time the space coordinates of the foot may 
be computed from equation (4) as 



Equation (49) is the equation of motion of the jth foot which is used when the foot is free. 


DIMENSIONLESS EQUATIONS 

The purpose of this part of the paper is to convert the equations of motion and cer- 
tain auxiliary relations into equivalent dimensionless forms. Working with dimension- 
less equations of motion facilitates the application of results to both model and full-scale 
versions of a vehicle as is frequently necessary in studies of lunar landing dynamics. 
Also, replacing time by a dimensionless variable allows one to rely to some extent on 
previous experience in sizing the time step for numerical integration of the equations of 
motion. 


Definitions of Dimensionless Quantities 

The symbol for a dimensionless quantity is generally formed by adding the letter 
B to the subscript in a corresponding dimensional quantity. (The only exceptions to 
this rule are the dimensionless quantities defined by equations (56a) to (57) and quantities 
such as Uxj^jUyj ? k>^Z,j,k anc * ^X,j’^Y,j’Wz j which are dimensionless as 
defined.) 
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Time .- Time t is replaced by the dimensionless variable tg where 



and L represents some characteristic length. An asterisk denotes differentiation 
with respect to tg. 


( 50 ) 


Lengths, forces, and velocities .- Lengths, forces, and linear velocities are ren- 
dered dimensionless by dividing them, respectively, by L, Mg, and \|gL. For example, 


Lengths: 

Forces: 

Velocities: 


XqB " L X 0 


F - F 

XB M g X 

* 1 
X™ . = — — x„ . 
FB ,] F ’J 


(51a) 

(51b) 

(51c) 


The dimensionless quantities arising in this manner (and denoted by subscript B) 
are as follows: 


From lengths 

From forces 

From velocities 

x B> y B’ Z B 

f xb» f yb> f zb 

V OXB’ V OYB’ V OZB 



* * * 

^B ,?, B’^B 

F HXB, j,k’ F HYB,j,k’ F HZB, j,k 

XOB’^OB’ Z OB 

X OB’ Y OB’ Z OB 

F HB,j,k 

S B,j,k 

* * * 

^PB’^PB’^PB 

F RB,j,k 

X HB,i,k> Y HB,i,k- Z HB,j,k 

X PB> Y PB’ Z PB 

F SB,],k 

S ECB,j,k 

* 

X FB,i> Y FB,j’ Z FB,j 

F ECB,j,k 

®EEB,j,k 

X HB,j,k' Y HB,j,k’ Z HB,j,k 

F EEB, j,k 

S CCB,j,k 

* 

*HB,j > k’%B > j,k’ ? HB,j,k 

F CCB,j,k 

S CEB,j,k 
* * * 

S B,j,k 

F CEB,j,k 

X FB,j’ Y FB,j’ Z FB,j 

S EB,j,k 

F H|B,j,k' F H77B,j,k’ F H?B,j,k 


S RB,j,k 

F FNSB,j 


S OB,j,k 

F FLXB,j> F FLYB,j’ F FLZB,j 


s SB,j,k 

F FLNXB,j’ F FLNYB,j’ F FLNZB,j 


W B, j 
H B,j 
D B,j 

F FLTXB,j ,F FLTYB,i ,F FLTZB,i 


^FB,j’^FB,j’^FB,j 
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Additional quantities.- Miscellaneous additional dimensionless quantities necessary 
for remaining developments are defined by the following equations: 


^ ^ ^ \ / * • *\ 
V OXB> V OYB» V OZB = gt V OX’ V OY> V Ozj 

(52) 

^B’^B’^B = \H( W £’ W 77 

(53) 

W £B’%B’^B = 

(54) 

T /N. N„ N-\ 
N N N - i± Is V S\ 

w «B»V ,N CB-g^ipyl^j 

(55) 

I 77^B’ I ?|B = ^(V I ?) 

(56a) 

%B^77B = ^(V^) 

(56b) 

^CB’Vb 

(56c) 

B ^B’ B r 7 B’ B ?B- ML 2( I rV I ?) 

(57) 


i /" 1 

P E,Z,j,k - Mg P ^,j,k ^ = 1 » 2 > 3 > 4 ) ( 58 ) 


K, 


SB 


L K 

Mg S 


L ' K i i 
K = ^ 

B »i»i Mg 


R NB,j> R TB,j 

E KB 



(59) 


(i = 1, 2, 3) (60) 


(61) 

(62) 
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Dimensionless Form of Equations of Motion of Body 


The following dimensionless form has been adopted for the equations of motion of 
the body. These dimensionless equations are readily derived from equations (9a) to (12) 
with use of the definitions of dimensionless quantities which have been established in 
this part of the paper. 


* 

XOB = V OXB 
★ 

y OB = V OYB 
* 

z OB = V OZB 
* 

V OXB = F XB 
* 

V OYB = F YB 
* 

V OZB = F ZB ' 1 

W £,B~ ‘‘SjB^B^rjijB " X ^b) + N ^B 
W t?B = w ?B W |B^77B ‘ ^|tjb) + N t)B 
^B = "SB^B^B ' Vb) + n cb 


<f> 
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cos 0 cos xjy 
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XPB " 

*ob 
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^PB 

Y PB " 

y ob 

)- 

6 p,q 
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^PB > 

Z PB " 

Z OB 
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J 

_ ..j 


V J 


(63a) 

(63b) 

(63c) 

(63d) 

(63e) 

(63f) 

(64 a) 
(64b) 
(64c) 


(65) 


( 66 ) 


Dimensionless Form of Equations of Motion of Feet 

The dimensionless forms adopted for the equations of motion of a penetrating foot 
and a free foot are, respectively, 
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X FB,j’ Y FB,j’^FB,j = R NB,j F FNSB,j( W X,j’ W Y,j’ W Z,j) 

+ R NB,j ( F FLNXB,j’ F FLNYB,j’ F FLNZB,j) 

+ R TB,j ( F FLTXBJ’ F FLTYB,j’ F FLTZB,j) ( 67 ) 


f ' \ 


- 

-i 

T 

f>© 1 

^FB,j 


r ~\ 

X FB,j 






XOB 

<Y FB,j 

>- 

6 p,q 


< 

^FB,j 

> + < 

y ob 

ZFB ,jj 





c (f) 

l ?FB »jJ 


Z OB 

^ J 


These equations are derived from equations (47) and (49). 

Dimensionless Form of Auxiliary Relations 


To facilitate working with the equations of motion in dimensionless form, a number 
of the relations established in the preceding sections are rewritten here in equivalent 
dimensionless forms: 

From equation (13): 


E KB 2 ( V OXB + V OYB + V SzbK( B ^!b + B t7B“^B + B £B w $b) ( 69) 


From equation (14): 




( X HB, j,k " X FB,j) 2 + ( Y HB, j,k " Y FB,j) 2 + ( Z HB,j,k " Z FB,; 


1/2 


(70) 


From equation (15): 


B,j,k s 


B,j,k 


X HB, j,k ■ X FB,j)( X HB,j,k " X FB,j 


+ ( Y HB,j,k " Y FB,j)( Y HB,j,k " Y FB,j) + ( Z HB,j,k " Z FB,j)( ! 


’ ^T7'T-» * I 1 Z TX n * 1. “ Z 


k FB,jll HB,j,k ” FB,j 


(71) 
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From equation (16): 


U X,j,k’ U Y, j,k’ U Z,j,k - sT~T ( X HB,j,k “ X FB,j)’ ( Y HB,j,k ‘ Y FB,j) ( Z HB,j,k " Z FB, j 

D 9 J J K L 


From equation (17): 


F HXB,j,k’ F HYB,j,k’ F HZB,j,k " F HBJ,k( U X,j,k’ U Y,;j,k> U Z,j,k) 


From equation (18): 


From equation (19): 


From equation (22): 


S SB,j,k - S RB,j,k " S B,j,k 


F HB,j,k " F SB,j,k + F RB,j,k 


F SB,j,k = P B,l,i,k + P B,2 


,j,k( S RB,j,k “ S B,j,k) + P B,3,j,k( S RB,i,k “ S B,j,k 


r 

+ P B,4,j,k ( S RB,j,k " S B,j,k)^ 


0 = S RB,j,k _ S B,j,k = S EB, j,k’ 0 < S B,j,k) (76) 


From equation (23): 


F SB,j,k - ” K SB ( S B, j,k ' S OB,j,k 


(° S S B,j,k ' S OB,j,k; ° <S B,j,k) (77) 


From equation (24): 


j=4 k=3 

f f f = y y ( 

r XB’ r YB’ r ZB L L\ 
3=1 k=l 


F HXB,j,k’ F HYB,j,k’ F HZB,j,k 


From equation (25): 


H£B,j,k 


HXB,j,k 


HT?B,j,k ; 


HYB,j,k/ 


H^B,j,k 


HZB,j,k 


32 



From equation (26a): 


j=4 k=3 

N ^ B = X X B7^( 1? HB,j,k F H?B,j,k ‘ ^HB,j,k F H?jB,j,k) 
j=l k=l s ' 

From equation (26b): 

3-4 k=3 

N '/B = X X B^( C HB,j,k F H4B,j,k - ^HB,j,k F H?B,j,k) 
j=l k=l ' i D ' 

From equation (26c): 

j=4 k=3 



From equation (28): 


W X,j X B + 


W YJ Y B 


W Z,j Z B + 


W B,j - 0 


From equation (30): 


H 


B 


- W X,j X B 


+ W Y,j Y B + 


W Z,j Z B 


+ W 




From equation (38): 


F FNSB,j 


K r i . 

Bj 1)3 B) j 


+ K n 9 D d 

B,2,3 B,j 


+ K B,3,j D 


3 

B,j 


From equation (39): 


D 


B,j 


W X,3 X FB,j 


W Y 
Y,i Y FB,j 


W Z,3 Z FB,3 


+ w 


B,j 


From equation (44): 


k=3 

F FLXBJ’ F FLYB,j’ F FLZB,j = ’ X F HXBJ,k’ F HYB,j,k’ F HZB,j,k 


(80a) 


(80b) 


(80c) 


(81) 


(82) 


(83) 


(84) 


(85) 
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From equation (45a): 


r 




r > 

F FLNXB,j 




F FLXB,j 

F FLNYB,j 

> = 

N 

p,q,3 

) 

\ 

F FLYB,j ? 

F FLNZB,j v 


L_ _j 


Jflzb,]^ 


From equation (45b): 


c 




r > 

F FLTXB,j 




F FLXB,j 

F FLTYB,j 

)- 

T p,q,i 


F FLYB,j / 

F FLTZB,i 




F FLZB,i 


(86a) 


(86b) 


NUMERICAL INTEGRATION OF EQUATIONS OF MOTION 

A brief discussion of the basic ideas involved in the numerical integration of the 
equations of motion is given in this section. 

Recurrence Equations 

The fundamental equations used in the integration are the following equations 
which form a set of recurrence equations. A superscript (n) indicates that a quantity 
is computed by or is otherwise associated with the nth of successive applications of the 
recurrence equations. The quantity At^ is the increment of dimensionless time used 
on the nth application. 


v-(n) Y (n) 7 (n) _ L(n-1) Y (n-1) „(n-l)\ (n)/*(n-l) *(n-l) *(n-l) 

*OB’ Y OB ,Z, OB " l A OB ’ *OB ,Zl OB / + Al B l A OB ’ OB ’ OB 


v (n) v (n) y (n) _ (y(n-l) y (n-l) y (n-l)\ t (n)/*(n-l) *(n-l) *(n-l) 

V OXB’ V OYB’ V OZB " \ V OXB’ V OYB’ V OZBy + Al B V v OXB ’ v OYB ’ v OZB 

Jn) (n) (n) _ / (n-1) (n-1) >-l)\ t (n)/*(n-l) *(n-l) *(n-l)' 

^B’^B’ £B “ ,w rjB ’ W £B y +At B ^ £B ’ rjB ’ W £B 


t( n ) vM jz( n ) - M n ) y( n ) y( n ) 
^B’* OB’ OB “ V OXB’ v OYB’ v OZB 


y( n ) y( n ) y( n ) - ( p( n ) p( n ) p( n ) _ 

V OXB’ v OYB’ v OZB - \ XB’ r YB’ r ZB / 


(87 a) 

(87b) 

(87c) 

(88a) 

(88b) 
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(89a) 


* (n) = (n) (n) /j \ N (n) 


*(n) _ w (n) (n) /j \ . (n) 

"sb^b^b 1 ^b; + N 7 jb 

*(n) = (n) (n)/j _ j \ N (n) 

?B W £B t/B(^£B Vb) ^B 


If foot j is penetrating at t^ 1 


x (n) y (n) z (n) = / x (n-l) y (n-l) z (n-l)\ ,(n)/* (n-1) * (n-1) * (n- 1) 

A FB,j’ Y FB,;j’ Z FB,j V X FB,j’ Y FB,j» Z FB,j; + At B ^ X FB,j’ Y FB,j > Z FB,j ; 


r( n ) Y"( n ) 7( n ) _ p p( n- l) Air W W \ 

FB,j’ Y FB,j’ Z FB,j - R NB,j F FNSB,j ( W X,j’ W Y,j’ W Z,jJ 


(n-1) 


r (F^ 11- ^) p( n- l) p( n 1) 
i'JBJ l r FLNXB,j’ r FLNYB,j ,r FLNZB,j 


+ • ( f{?- - —) F( n_1 ) fO 1-1 ) 


F' ' F' v pv 11 - 1 / 
TB,j \ FLTXB.j’ FLTYB,j’ FLTZB,j 


If foot j is free at tg 1 ^ 


[>) 1 
FB,j 


<Y (n) 

VFB,j 

■ 

z (n) 

l Z FB,t 

- 

ft(n) 1 
X FB,j 

|- 

( Y (n) 

Vfb,] 

II 

£( n ) 

l Z FB,jJ 



5(n) 

PA 


T r i 

t(n-l) 
^FB,j 


Wfb!!) * 


X OB 


r(n) 


f>(n-l) 
Ofb,^ 


An) 

J OB 


*(n) 

°p,q 


A T 

r ^ 

t(n-l) 

^FB,j 


r "\ 

y( n ) 

v OXB 

< 

Jn-1) 

^FB,j 

> + < 

y( n ) \ 

v oyb/ 


r(n-l) 


y( n ) 


I^fbJ 


t V OZBj 


If foot j is free at tj^ and free at t^ 1 "^ 

t(n) (n) Jn) _ ,(n-l) (n-1) r (n-l) 

^FBJ^FBJ^FBJ " ^FB,j^FB,j^FB,j 


(89b) 

(89c) 

(90) 

(91) 

(92a) 

(92b) 

( 93 ) 
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If foot j is free at t^ and penetrating at t^ 1 " 1 ^ or if foot j is penetrating at t^ 


ft(n) 1 

^FB,j 




f x (n) 

X FB,j 

x (n) 

■^OB 

^FB,j 

> = 

6 (n) 

P,q 

J 

Y ( n ) 

FB,j 

- Y (n) 
OB/ 

p (n) 

CFBjJ 




z (n) 

- 




L FB,j 

^obJ 


Auxiliary equations, for example, those relating forces and torques to the system 
variables, are used in the integration. Also, on each step a series of tests is performed 
to establish the dynamic state of each strut so that operations such as selection of appro- 
priate force equations and proper setting of reference lengths may be performed. The 
auxiliary computations and tests are not shown here but are described in complete 
detail in the section "Programing Instructions." 

Relation to Euler’s Method 

The integration scheme is practically identical to the classical straightforward 
method of Euler. (See ref. 9.) It differs from Euler’s method only in the following 
respect: 

Euler’s recurrence scheme, stated for simplicity for a single-degree-of-freedom 
system with dependent variable a 

“ = f( a »t B ) (95) 

amounts to the following: 

>) = a (n-l) + At W*(n-l) (96a) 

# n > = f(c/ n >4 n) ) (96b) 

However, if rate-dependent forces are considered, the single-degree-of-freedom analog 
of the dimensionless equations of motion does not take the form of equation (95) but is 
rather 

a= f(o,o,t B ) (97) 

* 

The variable a on the right appears because the forces and torques may depend upon 

the rates of change of the strut lengths. Thus, Euler's method is not directly applicable 
* 

because a is not an explicit function of a and tg. 
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To get around this difficulty and yet preserve the basic simplicity of Euler’s 
method, the dimensionless foot velocities . for a penetrating foot 

are not computed as functions of the forces on the foot at tj^ as would be required by 
a strict application of equation (96b) but as functions of the forces at t^ _1) as equa- 
tion (91) shows. Thus, the dimensionless accelerations and 

OXB’ OYB’ OZB 

^B’^B’^B are P arti ally dependent upon conditions at t^ 1-1 ^ because the forces 

F XB’ F YB ,F ZB and the torc l ues N |b ,N Sb ,N [b de P end on j’ ^FB j’ and ^FB j‘ 
In contrast to Euler’s method therefore, the recurrence scheme depends on the informa 
tionatboth t^ 1 and tg in computing the system variables (x^1,y1^,Z^), 

(y( n ) y(n) v (n) \ d / (n) (n) (n)\ 

\ OXB’ V OYB’ V OZB/’ and ^|B’ jjB’^CB/- 

If rate-dependent forces were suppressed, that is, if 


' EC,j,k’ F EE,j,k’ F CC,j,k ,F CE,j,k “ °’ 0 ’ 0 ’ 0 


(98) 


and if on the right in equation (91) the superscripts n-1 were replaced by n, the 
recurrence scheme would reduce to Euler’s method. 


Efficiency and Validity 

The method for integration was selected because it is very easy to program and 
because all the mathematical relations in the scheme may be interpreted physically. 
Efficiency was not a primary consideration. However, discussions with other workers 
in analysis of lunar landings have led the authors to believe that the method is relatively 
efficient in regard to consumption of computer time compared with methods now in use. 
About 2 minutes are required to compute an impact history with an IBM 7094 digital 
computer. Although the recurrence equations as given indicate a variable time interval 

fa) 

Atg , the authors to date have programed a constant time interval throughout an impact. 
The computer time could probably be reduced by using a short time interval during the 
initial part of an impact when accelerations are generally high and a longer one later 
when accelerations have been reduced. 

The authors are unable to offer any rigorous demonstration of the validity or accu- 
racy of the integration method. The method has been tried on a number of systems 
selected so that exact solutions for the motion could be computed, falling bodies, for 
example. The approximate solutions have consistently been in good agreement with the 
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exact when sufficiently small increments of the dimensionless time are used. Further 
confidence in the method has developed as a result of the good correlation between ana- 
lytical and test results reported in reference 1. 

PROGRAMING INSTRUCTIONS 

The object in this part of the paper is to provide instructions which expedite pro- 
graming a digital computer to carry out the computation of an impact history. The 
instructions are based entirely on relations which have been established in preceding 
parts of the paper. A constant time interval is assumed. 

Basic Organization 

The suggested basic organization of computing is shown in block diagram form: 



A table of input information is read which describes the vehicle, the landing sur- 
face, and the initial orientation and velocities of the vehicle. The input is converted to 
dimensionless form by use of the definitions and relations in section "Dimensionless 
Equations." Then initial values of all quantities of interest are computed. After ini- 
tialization, the computation proceeds into the integration cycle. On each pass through 
the integration cycle, the recurrence equations discussed in "Numerical Integration of 
the Equations of Motion" are applied to advance computation of the variables of the 
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system one step in time. As often as desired, quantities are extracted from the integra- 
tion cycle as output. Sometimes, it is convenient to convert the output to dimensional 
form. 


Input 

The quantities required as input are listed. The subscripts range as follows: 
3 = 1,2, 3, 4; k = 1,2,3; 1 = 1,2, 3, 4; i = 1,2,3; p = 1,2,3; q = 1,2,3. 

Gravitational constant: g 

Characteristic length: L 

Mass: M 


Principal inertias: 

Body coordinates of hard points: j )k ,< H>; j >k 

Initial body coordinates of feet: £^.,77^., 

* >3 f »] * >3 

Shock absorber constants: 

S E,j,k 

F EC,j,k ,F EE,j,k’ F CC,j,k ,F CE,j,k 
S EC,j,k’ S EE,j,k’ S CC,j,k’ S CE,j,k 


K S 


Boundary plane coefficients: A 
Surface impedance coefficients: 


X,j ,A Y,j’ A Z,j» A j 

R N j, R-p j, and Kjj ^Do not set R^ j or 

R-p j precisely equal to zero. In the present 
computing scheme, this procedure will cause 
the feet to be locked onto the landing surface 
plane. If infinite viscosity is desired, use very 
small finite values for R^ , and R-p 


Initial space coordinates of center of gravity: 



Initial space components of velocity center of gravity: 


v (°) v (°) v<°) 
v OX ,v OY’ v OZ 


Initial direction cosines of body axes with respect to space axes: 



39 



Initial components of angular velocity referred to body axes: 

Dimensionless time increment: Atg 

Conversion of Input to Dimensionless Form 
The computations necessary to render the input dimensionless are as follows: 






^B’Vb = 


^HB,j,k^HB,j,k’ ? HB,j,k " L (%,j,k’ 7? H,i,k^H,j,k) 


f (0) n (0) JO) _l(MJ0)J0) 

^FB,i’ r7 FB,i^FB,j " L^F^FJ^F,] 


l-l 


P, 


B,Z,],k Mg ^i, k 


S EB,j,k - I S E,j,k 


F ECB,i,k’ F EEB,i,k’ F CCB,j,k ,F CEB,i,k Mg ( F EC,j,k’ F EE,i,k» F CC,i,k’ F CE,j,k 


* * * * 1 /. * * • \ 
S ECB,3,k’ S EEB,i,k’ S CCB,j,k’ S CEB,i,k = ( s EC,j,k' S EE,i,k’ S CC,j,k’ S CE,i,k 


k sb = 

— K s 
Mg b 


R NB,j : 

,R TB,j = 

Ml If- 

¥ h 


L*K- . 


Kp -j -i 

Bjij] 

= hi 

Mg 


X (0) y(0) z (0) 

a OB’ 1 OB’ OB 

■i( 


y(0) y(0) y(0) _ * (y(0) y(0) y(^) 

v OXB’ v OYB’ v OZB “ ^\ v OX’ OY ,v OZ 
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JO) JO) (0) _ 

£B’%B’ ?B “ 





Integration Cycle 

The integration cycle is discussed before initialization because the necessity for 
the steps in initialization is much easier to understand once the integration cycle is 
understood. Upon entering the integration cycle for the first time, all quantities with 
superscript zero will have been provided as input or computed under initialization. 
Also the following quantities will have been computed under initialization: 


W X,j> W Y, i > W Z : 


r~~ — l 

N p,q,j 

and 



Upon entering the integration cycle for the nth time (n > 1), all quantities with superscript 
(n-1) will have been computed on the previous pass through the integration cycle. 

The organization of the integration cycle is shown in the block diagram which fol- 
lows. The circled numbers indicate the connections to be made from page to page. 

After the block diagram are listed the relations and tests necessary to perform the 
operations called for in each block. A block is identified by the number at the upper 
left-hand corner. 
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2 


For kth strut jth leg compute 
space coordinates of 
hard point 


For kth strut jth leg compute 
space components of velocity 
of hard point 


k = 1,2,3 


(12) i 

Test: Was jth foot free or 
penetrating at 
t (n-l) ? 


Penetrating - 


-Free 


Compute space coordinates of 
jth foot using equations for 
penetrating foot 


Compute space components of 
velocity of jth foot using 
equations for penetrating foot 


For kth strut jth leg compute 
strut length 


k - 1,2,3, 


Compute space coordinates of 
jth foot using equations for 
free foot 


Compute space components of 
velocity of jth foot using 
equations for free foot 

(27) r ~~~ 

For kth strut jth leg set 
strut length equal to its 

value at 


© 
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0 


j = 1,2, 3, 4 


3 = 1 , 2 , 3,4 


44 












0 


3 = 1 , 2 , 3, 4 
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BLOCK 1: 


f( n ) _ t (n-l) 

L T 3 — 


B 


+ At 


B 


BLOCK 2: 


X&) Yfe) zL n l =fx (n_1) Y (n '^ fvC 0 - 1 ) vfa" 1 ) v(n-l) 

vJB OB’ OB \ OB ,Y OB ,Z OB J + ^bI^OXB’^OYB ’^OZB 


BLOCK 3: 


v£0 y(n) y(n) _ fy(n-l) (n-1) v (n-l)\ ft-(n-l) *(n-l) *(n-l)' 

OXB’ V OYB’ V OZB \ OXB ’ V OYB ’^OZB/ + At B \ V OXB ’ V OYB ’ V OZB 


BLOCK 4: 


co 


^B^B - HbMb ’ ‘’.-fe 1 ’) ^ “b(^ 1) ^ 1) .^- 1) ; 


BLOCK 5: 


/■V-V") = (*<■>- V"- 1 V n - 1 >) + a , B (;<"- i),;(„-d) 


BLOCK 6: 




sin \p^ 


sin 0( n ) 

0 (n) 

> = 

cos \p^ 

$< n ) 


-cos 0^sin 


sin $( n ) 


cos 1 
sin $( n ) 

0 


-8 

-sin i// n ) 

0 

( 

J n V 



VB 

cos 0^cos \ // n ) 

1 


JK) 

W ?B 

L. J 

sin ^( n ) 



BLOCK 7: 


E 


(n) _1 
KB "2 







+ B 


rjB 



+ b^ b 



BLOCK 8 : 

6 1^1 = cos ^^ cos 0^ - cos 0^sin (/>^sin i// n ) 
2 “ sin 0^sin 0^ 
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BLOCK 8 - Concluded: 


6^ n l = -cos p^ sin p^ - cos fi^cos p^sin p^ 
1,3 

= - sin p^cos p^ - cos 0^sin ciscos p^ 


>> 

5 2,2 

= COS 

xp^ sin 

e ( n ) 

.(n) 

6 2,3 

= sin 

xp^sin 

pW 

c( n ) 

6 3,1 

= sin 

0^sin 

p( n ) 

,(n) 

6 3,2 

= cos 

0 (n) 


*(n) 

6q 0 

= sin 

0 {n) cos 

p<“) 


BLOCK 9: 


tin) 

1,1 


cos p^sin 4>i n ^ + cos f/ n ^cos P^sin p^ 


;<«> 


+ I sin O^sin p^sin p^ 


?w 


S< i n, 2= + 


- sin p^cos (pi n ^ + cos 0^sin P^cos 1 pi n ^ xj/i n ^ 

L J 

sin p^cos 0^1 0^ + 


. (n) . „(n) [ * (n) 

cos p' 7 sin 0 ' '\ p' ' 


tin) = 


1,3 


cos p^cos <pi n ^ - cos 0^sin p^ sin p^ n \ p^ 


sin 0^cos p^sin p^j0^ 


+ sin p^sin p^ - cos 0^cos p^ cos p^ n i 


$< n > 


6^ = isin p^sin p^ - cos 0^cos p^cos p^ 

2.1 l T • 


}W 


+ [sin (/ n ^sin p^cos p^ n i 0^ 

- cosp^cos p^ - cos 0^sin p^sin p^ 


pW 
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BLOCK 9 - Concluded: 


BLOCK 10: 


BLOCK 11; 


S^ 2 = [cos 0^cos _ sin ^( n ) sin e ( n ) 




= _ sin ^ (n) cos 0^ + cos 6>^sin 0^cos 0^ 

+ |sin e^cos ^Wcos 0 (n ^ (n) 

+ [cos 0^sin 0^ + cos 0^cos 0^sin 0^ 






6g n j = [sin fl^cos 0^ 


0^ + cos (J^sin 0^ 


S( n ) 


5 *(n) _ 


3,2 


r 

[sin 6 


(n) 


5(n) 


t(n) - 


3,3 


sin e^sin 0^ 


0 (n) + 


cos 0^COS 0^ 


§(n) 


r(n) 
HB,j,k 


,(n) 

J HB,j,k 


t 

r 

' 


T 

r > 

*HB,j,k 


r x ( n )] 

X OB 

j- 


6 (n) 

p,q 

< 

j^HBJ.k 


Y ( n ) \ 
*ob/ 

w 

J 




jHB,j,k 


7 (n) 

Z OB 


fu n) 

^HB,j,k 


Jy^ 


\ HB,j,k, 

) 

z (n) 


1 HB, j,k 



S( fi ) 

p,q 


r 


HB,j,k 
{ ^HB ,j,k 
C HB,j,kj 


+ 


v (n) ^ 

OXB 

V (n) 

V OYB/ 

y( n ) 

v OZB 

v j 


BLOCK 12: The foot was free if either H^" 1 ) i? 0 or V^" 1 ) > 0; otherwise, the foot 
was penetrating. ’ ,J 
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BLOCK 13: 


(n) (n) (n) 

A FB,j’ X FB,j’ rB,j 



Y ( n ' !) 
,Y FB,j 


,Z 


(n- 1)\ 
FB,j J 


+ At B 


/*(n-l) * (n-1) * (n- 1)\ 
\ X FB,j’ Y FB,j^FB,jj 


BLOCK 14: 


*(n) * (n) i(n) 

A FB,i’ 1 FB,j’ r B,j 


R NB, i F ™SB,i( W X,i’ W Y, i ’ W Z,i) 

+ R NB t )( F FLNXB,]’ F FLNYB,j ,F 

, r f F ( n_1) F (n_1) F 

+ K TB,jr FLTXB,j’ r FLTYB,i’ r 


(n-1) \ 

FLNZB,jy 

(n-1) \ 

FLTZB, jj 


BLOCK 15: 


S 


(n) 

B,j,k 


(n) „(n) 

■HB,j,k ' A FB,j 


.(n) 

HB,j,k 



4 - 


iz (n) 

\ Z HB,j,k 



1/2 


BLOCK 16: 





(n) 

HB,j,k 



*(n) 

HB,j,k 



BLOCK 17 : 

Jn) _ (n- 1) (n) 

b SB,j,k " b RB, j,k ' S B,j,k 


BLOCK 18: The following table gives conditions and the operations to be performed when 

all conditions in a row of the table hold. Test first to determine whether the first condi- 
tion holds; all other conditions are mutually exclusive. 
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Contracted 

4t, k >o 


Contracting I 

S^ n) . <0 

B,j,k 


Conditions 


Elastic 


0 - c( n ) - o 

& SB,j,k b EB,j,k 


o • s (n) - s (0) 


Do operation 


„(n) _*(n-l) 

°RB,j,k “ b RB,j,k 


r (n) _ K f JO) j n ) \ 

SB,j,k “ ^SB^ S B,j t k ' S B,j,kj 


S n , ^ <S, 


F fn) - 0 
* RB,j,k ' U 


cW _ c(n- 1) 
a RB,j,k " b RB,j,k 


M 


B,j,k| CCB,j,k I r SB,j,k : 


F <“) - . 
RB .j»k “ J 


S(n) 

B,j,k 

a CCB,j,k 


CCB,j,k 


,(n) <,(n-l) 

s RB,j,k RB,j,k 


Jr\) 

~B,j,k| ^CCB.j.k r SB,j,k = 


| S B,j,k| < S CCB,j,k 


>) F 

RB,j,k ” r CCB,j,k 


<>) Jn-l) 

S RB,j,k * b RB,j,k 


F (n) -0 
* SB,j,k “ 0 


» 


*(n) 

^B,j,k 

S CCB,j,k 


CCB,j,k 


Yes 

No 

1 1+ \ + 

- S 

1 

1 

1 | B,j,k| & CCBj,k 


s (n) _ (n-1) 
kB, j,k ~ b RB,j,k 


F (n) - 0 
*SB,j,k " U 


(n) 

F RB,j,k = F CCB,j,k 


„(n) _ (n-1) 

RB,j,k = nB,j,k 


B,j,k ^CEBjj^k 




F SB,j,k - P B,l,j,k + P B,2,j,k( S SB,j,k) + P B,3,j,k( S SB, j,k) + P B,4, j t k( S SB,j,l 


p(n) _ _ 

r RB,j,k " T 


q(n) 


'CEB, j,k 


r CEB, j,k 
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Conditions 

Contracted I Contracting j Elastic 

si "2 - > 0 | ,, < 0 * 0 • 


SB,j,k 


B t j,k 


SB,j,k EB.j,k 


Do operation 


(n) Jn- 1 ) 

s RB,j,k “ RB.j.k 


s > s p(tt) 

B,j,k CEBJ,k r SB 

■ Jv) 


•- p +. j> { \ + p (s^ \ + p /g( n ) Y* 

,j,k " R,IM,k\' SB,j,kJ 1 B,3,j,k\^SB, j,kj + 1 B,4,j,k^SB,j,kJ 


' 1 


s (n) s (n) + s 

& RB, j,k b B,j,k + ft EB ( ].k 


F S»,j,k " P B,l,j,k + P B, 2 ,j,k( S EB,j,k) + P B, 3 ,j,k( S EB } 3 ,k ) 2 + P B, 4 ,j,k( S EB,j,k ) 3 


F (n) 

r RB,i,k ’ + 

S, 


g(n) 

B r j,k 


CEB, j, k 


! s (n) - s (,l) + s 

! RB,j,k ' B,j,k + EB,j.l 


S B,j,k ‘ S CEB, j,k F SB,j,k " P B.l,j.k + P B, 2 , j f J S EB,j,k) + P B, 3 t j,k( S EB,.i t kJ 2 + P B, 4 ,j.k( S EB,j,k) 

» 


F RB, j,k " " F CEB,j,k 


»(n) eO 1 ’ 1 ) 

| s RB,j,k HB,j,k 


B P j,k| ECB,j,k i 


r 0 

SB,j,k u 


F (n) - 
r RB, j,k 


*(n) 

b B,],k ] 
S ECB,i,k 


(n) „(n-l) 

s RB,.j,k RB,j,k 


No 


S B,j,k| * S ECB,j,k 


F (l ° ~ 0 


» 


B,j,k| < S EEB,j,k 


F RB,j f k F ECB,j,k 


«(n) (n- 1 ) 

RB,j,k RB,j,k 


F » = o 

SB,j,k u 


S B,j,k| ^EEB, j,k 


*(n) 

P r 
* RB, i.k + h EEB,i.k 

S EEB, i,k 


(n) ^(n-l) 

ft RB,i,k “ s RB,j,k 


F 0l) - 0 
SB,j,k “ u 


F 0l) -F 
r RB, j,k r EEB i,k 
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BLOCK 19: 


^j( n ) p( n ) TT^ - ^ 

U X,j,k’ U Y,j,k’ U Z,j,k (nj~ 1X 

S' . . C 
B ,j,k 


= p( n ) , F ( n ) 

HB,j,k SB,j,k + RB,j,k 


■ (n) - x<t ( Uy<"> 


-*£LU* W 


^(n) 

HB,j,k ^FBj/’^HBJ.k " *FB,jM A HB,j,k " Z FB, 


,(n) (n) (n) (n) / (n) (n) TT (n) \ 

HXB,j,k’^ HYB,j,k’ F HZB, j,k " F HB, j,k ^ U X,j,k’ U Y,j,k’ U Z,j,k) 


BLOCK 20: 


F (n) ^ 


“ 

— 


f F (n) 1 

H^B,j,k 





HXB,j,k 

J F ( n ) 

H7,B,j,k 

> = 

6 ( n ) 

p,q 


< 

F ( n ) 

HYB,j,k 

F (n) 

L HCB,j,kJ 





F (n) 

^HZB,j,kJ 


BLOCK 21: 


(n) (n) (n) 

HXBT,j’ r HYBT,j’ 1 HZBT,j 


k=3 


k=l 


p( n ) F 0^) p(^) 

HXB,j,k’*HYB,j,k’*HZB 


,j,k) 


BLOCK 22: 


f F ^ n ) ^ 

%B,j,k r H77B,j,kJ 


, v k=3 

N ( „ n) - Y 1 ( n >) 

|BT,j 4 B^ B rHB,j,k F H?B,j,k “ 

K — 1 
k -3 

N (n) . = t_i Jr >) _ t » \ 

??BT,j BjjB \^HB,j,k h£B,j,k %B,j,k F HCB,j,k j 


_ X' 1 ( i- -p. C n ) _(n) 

? BT >j - 4 ^HB, j,k F H7?B, j,k " ^ ^ 


* B B, j,k H£B, j 


,j,k) 


BLOCK 23: 


= 

B,j 


Kj)( x fL 


r(n) 

'Y.jA'FBJ 


+ (W v JlY^A J + [ W, 


:,j)( z FB,j) 


+ W B,i 


BLOCK 23 - Concluded: 


F (n) 

FNSB,j 

BLOCK 24 : 

F ( n ) 

FLXB,j 


= K r 


+ K t 


d<">V 


+ K r 


B, l,j B, j + ~B,2,j\ B,jj T B,3, B,j 


D 


(n) 


(n) (n) _ / (n) (n) 

,r FLYB,j’ r FLZB,j “ ” y r HXBT,j ,r HYBT,j’ 


fF (n) 1 

* FLNXB,j 





r F (n) ] 

FLXB,j 

F (n) 

y FLNYB,j 

>- 


N p,q,j 

< 

f(°) V 

FLYB,j / 

(n) 

JfLNZB.jJ 





F (n) 

L FLZB,j J 


r F (n) i 

r FLTXBJ 




f F (n) 1 

FLXB,j 

J F (n) 

Y FLTYB.j 

> = 

T p,q,i 

< 

F (n) 

r FLYB,j / 

F (n) 

FLTZB,j 



1 

F (n) 

If FLZB,j J 


BLOCK 25: 


o> 1 

X FB,j 



T 

ft n-lfi 
^FB,j 


x (°r 

c)B 

y( n ) 

Y FB,j 

= 

6 (n) 

p,q 

< 

>-!) 

^FB,j 

> + < 

Y (n) } 
*OB 

7 (n) 

Z FB,j v 




r (n-l) 

OFB,^ 


7 ( n ) 

L z obJ 


BLOCK 26: 


*(n) 

X FB,j 




T 

ftn-lfl 

^FB,j 


r (°) i 

V OXB 

iv (n) 

SY FB,j 

- 

*(n) 

6 p,q 


< 

(n-1) 

^FB,j 

> + < 

V (n) 

v OYB 

£(n) 

l Z FB,t 





>-1) 

L ? FB,t 


v (n) 

LozbJ 


,(n) 

HZBT, j 
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BLOCK 27: 


BLOCK 28: 


BLOCK 29: 


BLOCK 30: 


BLOCK 31: 


BLOCK 32: 


BLOCK 33: 


BLOCK 34: 


BLOCK 35: 


BLOCK 36: 


c(n) o(n-l) 
B,j,k - b B,j,k 



,(n) _ (n-1) 

SB,j,k " oB,j,k 


i( n ) _ c(n) 
RB,j,k - b B,j,k 


(n) (n) (n) 

HXBT,j ,r HYBT,j’ r HZBT,j U ’ U ’ U 


N |BT,j’ N T?BT,j’ N £BT,j = °’ 0 ’ 0 


F' ' = 0 

FNSB,j u 


p( n ) p( n ) p( n ) 

FLNXB,j ,r FLNYB,p FLNZB,j 

p( n ) p( n ) p( n ) 

FLTXB, j’ r FLTYB, j»* FLTZB,j 


0,0,0 

0,0,0 


H B, } j = W X,j X FB,j + W Y,j Y FB, j + W Z,j Z FB,j + W B,j 




y-( n ) 

FB,j 


+ W. 


Y,j 


&(n) 

FB,j 


+ W r 




S(n) 

FB,j 
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BLOCK 37: The foot is free if either H^. = 0 or . > 0; otherwise, the foot is 

penetrating. 

BLOCK 38: The foot was free if either = 0 or > 0; otherwise, the foot 

h,] Nt>,j 

was penetrating. 

BLOCK 39: 


an) ^ 

^FB,j 




an) 

x FB,j ' 

Y (n)^ 

A OB 

>> 
y FB,j 

> = 

6 (n) 

p,q 

< 

Y (n) 

Y FB,j ' 

Y (n) 
*OB , 

r (n) 

l C FB,jJ 




7 (n) 

l Z FB,j ' 

7 (n) 

^OB, 


BLOCK 40: 


,(n) (n) (n) (n-1) (n-1) (n-1) 

^FB,j ,r ^FB, j’^FB,j ^FB,j’ ? 'FB, j’^FB, j 


BLOCK 41: 


TT (n) F (n) F (n) -Vf F (n) F (n) F (n) 
*XB’ YB’ ZB " L V HXBT,j’ HYBT, j’ r HZBT,j 


j=l 


BLOCK 42: 


A( n ) y( n ) y( n ) j-( n ) p( n ) p( n ) _ J 

v OXB’ v OYB’ v OZB XB’ YB’ ZB 
BLOCK 43: 


}=4 


r(n) M (n) *» 


N ^B ,N r;B ,N ^B " ^^ N ?BT,j» N iBT,j» N CBT,^ 


BLOCK 44: 


*(n) (n) (n), \ . „(n) 

W £B " W 77B w Cb(Vb " ^B) + N £B 


* (n) (n) (n) 


w i 7 B = w 


?B w |Bfci?B " %b) + N rjB 


* (n) (n) (n) 


(n) 

7 B 

(n) 


“?B = “^B^rjB^^B ' Vb) + N CB 
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Initialization 


Before entering the integration cycle, it is first necessary to compute certain 
constants and to compute initial values of those time-dependent variables whose initial 
values are not given as input. As in the previous section, the computation will be 
described by a block diagram followed by the details of the steps in each block. When 
referring to equations written in the description of the integration cycle it is always 
assumed that the superscript (n) has the value zero. 
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© 


( 11 ) 



© 







© 


0 


( 16 ) 



= 1 . 2 , 3.4 


60 


5 


© 












BLOCK 1 : 

t (0) - 0 
l B u 

BLOCK 2: 


B £B’ B i7B’ B ?B - 


BLOCK 3 : Use block 7 of integration cycle. 


BLOCK 4: Compute 


* 


( 0 ) 


= arc cosine 


6 


( 0 ) 

3 ^ 


sfi ) 2 + U }«)Y 


3,3/ 


(° = 


<$> 


( 0 ) 


5 77 
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Set = - 4>^ if 6 ^ <0; otherwise, set 


*<»> , W 


Compute 


Compute 


0 ^) - arc cosine 


= arc cosine 




2 . (M \ 2 


0 < e ^ < 


0 S |^| g v 


Set 4^ = - if 6 j 2 < 0; otherwise, set 




If 5 ^) = 0, then <f>^ = arc cosine (1) =0. However, because of round-off 

3 > 1 

error, it is possible for a computer to generate a number slightly larger than unity for 

h 2 

/U°l and thus cause an error stop in most subroutines for computing the arc 
3,3/ II 3,3 ’ 

cosine. This possibility should be accounted for in programing the computation of the 
initial Eulerian angles. 

BLOCK 5 : Use block 6 of the integration cycle. 

BLOCK 6 : Use block 9 of the integration cycle. 

BLOCK 7 : Compute j>Wy j’^z j ° r 
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BLOCK 8: Compute 


N 


p,q,j 


where 


Compute 


[ T p,q>j] 


where 


N l,l,j W X,j 

N l,2,r N 2,l,j =W X,i W Y,i 

N 1 ,3,r N 3,l,J= W X,i W Z,i 

N,, f =wj. 

2,2,3 Y,j 

N 3,3 (j =W 2 Z,i 


T 1,1,J " W Y,j + W Z,j 

T u,i = T M,rVY,i 
T l,3,r T 3,l,i- W X,j W Z,i 
T 2,2,j = W X,j + W Z,j 
T 2,3,i = T 3,2,j ' - W Y,j W Z,i 

T 3,3,J ' W l,j * 


BLOCK 9: 


fx (0) 1 
X FB,j 


J Y «» 

rFBj 

\ 

/ = 

z (0) 

L Z FB,t 



6 (0) 

p,q 


t (0) ' 
’FB,j 


J „(°) 
\^FB,j 


CjM] 

*OB 


!>♦ ^ <0) 


JO) 

Ofbj 


OB 




z (0) 

L z obJ 


BLOCK 10 : Use block 35 of integration cycle. 
BLOCK 11 : Use block 10 of integration cycle. 
BLOCK 12 : Use block 11 of integration cycle. 



BLOCK 13 : Use block 15 of integration cycle. 

BLOCK 14: 


S 


( 0 ) 

SB,j,k 


= 0 


BLOCK 15 : 

c(°) _ o(0) 

b RB,j,k “ & B, j,k 

BLOCK 16: If £ 0, the foot is to the surface side of the plane; otherwise, the foot 

B, j 

is not to the surface side of the plane. 


BLOCK 17: 


t (0) £(°) £(°) - 0 0 0 

X FB,j’ Y FB,j’ Z FB,j °’°’ U 

Use block 16 of integration cycle. 

Use blocks 18 to 22 of the integration cycle and delete from block 18 the 
(n) 

RB,j,k‘ 

Use block 23 of integration cycle. 

Use block 24 of integration cycle. 

Use block 26 of integration cycle. 

Use block 28 of integration cycle. 

Use blocks 31 and 32 of integration cycle. 

Use block 33 of integration cycle. 

Use block 34 of integration cycle. 

Use block 36 of integration cycle. 

Use blocks 41 and 42 of integration cycle. 

Use blocks 43 and 44 of integration cycle. 


BLOCK 18 : 

BLOCK 19 : 
calculation of reference length S 
BLOCK 20 
BLOCK 21 
BLOCK 22 
BLOCK 23 
BLOCK 24 
BLOCK 25 
BLOCK 26 
BLOCK 27 
BLOCK 28 
BLOCK 29 


Comments on Output 

Printing or plotting quantities generated by a digital computer can involve sub- 
stantial computer time. Therefore, when the procedure described in this paper is pro- 
gramed, outputting should be made flexible so that the number of output quantities can be 
kept to a minimum consistent with the objectives of a particular investigation. It is 
usually not necessary to output after each execution of the integration cycle since the 
time step required for satisfactory integration is usually much smaller than that required 
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for plotting. Conversion from dimensionless to dimensional output is readily programed 
by using the inverse forms of the relations given in "Dimensionless Equations." 

RELATION BETWEEN PRESENT PAPER AND REFERENCE 1 

Refeience 1 gives the results obtained by using the procedure to compute landing 
stability boundaries, that is, to compute bounds on approach velocities and orientations 
within which a vehicle will not overturn and beyond which it will overturn. The calcula- 
tions were made for a 1/6-scale dynamic model of a lunar landing vehicle suitable for 
manned landing. The computed stability boundaries were compared with boundaries 
obtained by landing tests of the model. 

The computed stability boundaries presented in reference 1 fall into three cate- 
gories termed "elastic shocks," "inelastic shocks," and "rigid body." The first two 
categories were obtained by using an IBM 7094 digital computer programed to execute 
a procedure in all essentials equivalent to the procedure of this paper. The boundaries 
designated "rigid body” were obtained by using a procedure reported in reference 5 
which is based on the assumption that the entire vehicle moves as a rigid unit. The dif- 
ference between the boundaries designated "elastic shocks" and "inelastic shocks" lies 
in the representation of the shock absorber. In the first instance, there is assumed to 
be a spring in the shock absorber; in the second, the shock is assumed to provide no 
elastic restoring force. 

For precise cross-referencing, the numerical values of input to describe the model 
and the landing surface of reference 1 are given in table I on page 66. 

SOME NEEDED IMPROVEMENTS IN ANALYSIS 

It is extremely important to render analytical procedures much more efficient in 
regard to consumption of computer time in order to consider the multitude of landing 
situations which may be encountered. Design of vehicles and planning of missions would 
be greatly enhanced if the time for calculating an impact history could be reduced to the 
order of a hundredth of a minute as opposed to present times of 2 minutes or more with 
fast computers. The possibilities for improvement in this area are largely unexplored. 

A systematic method is needed for including the effects of overall elasticity of the 
vehicle. A way should be devised so that the data describing the elastic characteristics 
of the system can be obtained either by structural analysis or by feasible tests if the 
vehicle is available. 
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TABLE I.- INPUT FOR LANDING STABILITY CALCULATIONS REPORTED IN REFERENCE 1 
= 32.2 ft/sec 2 (9.8146 m/sec 2 ); L = 2.26 ft (0.68885 m); Atg = 0.000369 secj 

(a) Input for configuration 


lb-sec^ 


2.08 ft-lb-sec ^ 
2.08 ft-lb-sec 2 

1.098 ft-lb-sec 2 
0.895 ft 
-0.370 ft 
-1.47 ft 
1.4675 ft 
0 

-0.647 ft 
0.895 ft 
0.370 ft 
-1.47 ft 
0.370 ft 
0.895 ft 
-1.47 ft 


1.4675 ft 
-0.647 ft 
-0.370 ft 
0.895 ft 
-1.47 ft 
-0.895 ft 
0.370 ft 
-1.47 ft 
-1.4675 ft 


-0.647 ft 
-0.895 ft 
-0.370 ft 


0.2876 m-kg-sec 2 
0.2876 m-kg-sec 2 

0.1518 m-kg-sec 2 
0.273 m 
-0.113 m 
-0.448 m 
0.447 m 
0 

-0.197 m 
0.273 m 
0.113 m 
-0.448 m 
0.113 m 
0.273 m 
-0.448 m 
0 

0.447 m 
-0.197 m 
-0.113 m 
0.273 m 
-0.448 m 
-0.273 m 
0.113 m 
-0.448 m 
-0.447 m 
0 

-0.197 m 
-0.273 m 
-0.113 m 
-0.448 m 


-0.370 ft -0.113 m 

-0.895 ft -0.273 m 


-1.47 ft 
0 

-1.4675 ft 
-0.647 ft 
0.370 ft 
-0.895 ft 
-1.47 ft 


-0.448 m 
0 

-0.447 m 
-0.197 m 
0.113 m 
-0.273 m 
-0.448 m 


2.205 ft 0.672 m 


-2.26 ft -0.689 m 


2.205 ft 0.672 m 


-2.26 ft -0.689 m 


-2.205 ft -0.672 m 


-2.26 ft -0.689 m 


-2.205 ft -0.672 m 


-2.260 ft -0.689 m 




TABLE I.- INPUT FOR LANDING STABILITY CALCULATIONS REPORTED IN REFERENCE 1 - Concluded 


(b) Input for landing surface 


A X,3 

0 

0 

A Y.i 

0.17364818 

0.17364818 

> > 
^ N 

0.98480775 

0.98480775 

0 

0 

“nj 

10 ^ ft/lb-sec 

0.672 x 10 ** m/kg-sec 

R T,j 

10 H ft/lb- sec 

0.672 X 10 11 rr/kg-sec 

K i,j 

0 

0 


(c) Input for shock absorber 


For elastic shocks: 


P l,j,k 

0 

0 

P 2,j,l 

53,600 lb/ft 

79,765 kg/m 

P 2,j,2 

73,900 lb/ft 

109,975 kg/m 

P 2,j,3 

53,600 lb/ ft 

79,765 kg/m 

P 3,j,k 

0 

0 

P 4,j,k 

0 

0 

S E,j,l 

0.00239 ft 

0.000728 m 

S E,j,2 

0.003465 ft 

0.00106 m 

S E,j,3 

0.00239 ft 

0.000728 m 

F EC,j,k 

5 lb 

2.27 kg 

F EE,j,k 

5 lb 

2.27 kg 

F cc,j,i 

128 lb 

58.1 kg 

F CC,j,2 

256 lb 

116 kg 

F CC,j,3 

128 lb 

58.1 kg 

F CE,j,k 

5 lb 

2.27 kg 

S EC,j,k 

0.01 ft/sec 

0.003048 m/sec 

S EE,j,k 

0.01 ft/sec 

0.003048 m/sec 

S cc,j,k 

0.01 ft/sec 

0.003048 m/sec 

S CE,j,k 

0.01 ft/sec 

0.003048 m/sec 

*8 

73,900 lb/ft 

109,975 kg/m 


For inelastic shocks: 


i,j,k 

0 

0 

S E,j,l 

S E,j,2 

S E,j,3 

0.00203 ft 

0.000619 m 

0.00350 ft 

0.00107 m 

0.00203 ft 

0.000619 m 

F EC,j,k 

5 lb 

2.27 kg 

F EE,j,k 

5 lb 

2.27 kg 

F cc,j,i 

128 lb 

58.1 kg 

F CC,j,2 

256 lb 

116 kg 

F CC,j,3 

128 lb 

58.1 kg 

F CE,j,k 

5 lb 

2.27 kg 

S EC,j,k 

0.01 ft/sec 

0.003048 m/sec 

S EE,j,k 

0.01 ft/ sec 

0.003048 m/sec 

s cc,j,k 

0.01 ft/sec 

0.003048 m/sec 

S CE,j,k 

0.01 ft/sec 

0.003048 m/sec 

K S 

73,900 lb/ft 

109,975 kg/m 


67 




For the usual vehicle configuration, a rocket nozzle protrudes downward beneath 
the body when the vehicle is in a landing attitude. As the shock absorbers stroke, the 
nozzle may impinge on the landing surface and affect the course of the vehicle during 
the impact and possibly affect the stability against overturning. The subject of nozzle 
impingement requires study and documentation. 

CONCLUDING REMARKS 

This paper has given the development of a procedure for computing the motions 
during impact of a spacecraft with legs representative of presently conceived vehicles 
for lunar landing. 

Idealization of the vehicle and landing surface is discussed in a general way. The 
vehicle is treated as an arbitrary rigid body to which there are attached up to four legs, 
each leg consisting of three struts in an inverted tripod arrangement. The struts are 
connected to the body by universal joints and the junction point of the three struts at the 
foot of a leg is also a universal joint. There is a shock absorber in each strut. The 
individual struts may shorten or lengthen because of stroking of the shock absorbers 
but otherwise do not deform. Locations of the points where the struts attach to the body 
and the initial positions of the feet relative to the body may be arbitrarily chosen. The 
legs are considered to have no mass. The shock absorber in a strut is considered to 
produce forces directed along the axis of the strut. The magnitudes of the shock 
absorber forces are assumed to depend on the instantaneous length of the strut and the 
rate of change of this length. The boundary of the landing surface is represented by a 
set of arbitrarily oriented planes, one plane associated with each foot. If a foot is inter- 
acting with the surface material, a force is assumed to act on the foot. This interaction 
force is assumed to be a function of the position and velocity of the foot. 

A derivation of the differential equations which govern the motion of the rigid body 
part of the idealized vehicle is given. The equations consist basically of Newton's equa- 
tions of translation of the body and the Euler equations for rotation of the body. Forces 
and torques from the shock absorbers appear in the equations as specified variables. 

Idealization of the shock absorbers is discussed in detail. Emphasis is placed 
on practical considerations in representing the aluminum honeycomb shock absorbers 
which were utilized on a model vehicle reported in NASA TN D-4215. These considera- 
tions include representing constant force crushing characteristic of honeycomb, repre- 
senting the overall vehicle elasticity by springs in the shock absorbers, and representing 
gross system damping to prevent unrealistic bounding of the idealized system. Expres- 
sions are derived relating shock absorber forces and torques to the system variables, 
and thereby completing the equations of motion of the body. 
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The idealization of the landing surface is developed and equations of motion of the 
feet are derived. In representing forces generated by the interaction of the feet with the 
surface, two facts must be considered. First, knowledge of the properties of the lunar 
surface which would affect landing performance is as yet limited. Second, soil mechanics 
has not progressed to the point where one can predict with any confidence the forces on 
an arbitrary body impinging upon or passing through soil even under laboratory condi- 
tions. Therefore, cases where the feet move substantially through or along the surface 
are studied by assuming laws for the force on a foot. The definitions and associated 
equations allow one to construct a variety of laws for the force on a foot by specializa- 
tion of constants. The derivation of the equations of motion of the feet rests on the 
assumptions that a foot not interacting with the landing surface material moves to main- 
tain a static balance between the interaction force and the forces from the shock 
absorbers bearing on the foot. 

The equations of motion of the body and feet and necessary auxiliary relations are 
converted to equivalent dimensionless forms. Working with the dimensionless equa- 
tions facilitates application of results to both model and full-scale versions of a vehicle. 
Also, replacing time by a dimensionless variable allows one to rely to some extent on 
previous experience in sizing the time step for numerical integration of the equations 
of motion. 

The scheme for numerical integration of the equations of motion is given; the 
method used is a slight modification of Euler's method. 

Explicit instructions are given for programing a digital computer to compute a 
general impact history. 

Some improvement is needed in the analysis of landing vehicles with legs. Ana- 
lytical procedures should be made much more efficient in regard to consumption of 
computer time, a systematic method is needed for including overall elasticity of the 
vehicle, and the effects of impingement of the rocket nozzle on the landing surface 
require study and documentation. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., May 8, 1967, 

124-08-04-13-23. 
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